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Abstract 



Millisecond crystal relaxation has been used to explain anomalous decay in 
doped alkali halides. We attribute this slowness to Fermi-Pasta-Ulam solitons. 
Our model exhibits confinement of mechanical energy released by excitation. 
Extending the model to long times is justified by its relation to solitons, 
excitations previously proposed to occur in alkali halides. Soliton damping 
and observation are also discussed. 
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Crystals respond in picoseconds. This is the characteristic time for all sorts of lattice 
adjustments, for example passage to a "relaxed excited state." Nevertheless, an anomaly in 
the decay of luminescence in doped alkali halides has led us to conclude that relaxation 
can sometimes take place on a scale of milliseconds. It's as if the tides moved like tectonic 
plates. In this article we propose that the formation and slow decay of Fermi-Pasta-Ulam 
(FPU) solitons |3|-[5| can account for this dramatic slowdown. The presence of FPU solitons 
in alkali halides has been previously proposed [EL|7| and we believe that particular features 
of our system |]|| enable their production and observation. 

We begin with an overview. KBr (say) is dilutely doped with Pb 2+ . One can think of 
the Pb and its 6 surrounding Br's as a quasimolecule. This quasimolecule is excited by a 
flash of UV light and because of the Jahn- Teller effect distorts considerably. In the absence 
of the constraining lattice one would expect expansion along one of the axes of perhaps 
15%. The distortion is asymmetric (there is shrinking along other axes) and the Jahn- 
Teller interaction causes well-understood degeneracy-breaking effects. It is the subsequent 
return to the ground state that provides evidence for slow lattice relaxation. The excited 
quasimolecule can go to a radiative level with a 25 ns lifetime or to a metastable level with 
an 8 ms lifetime. The luminescence decay however is not the sum of two exponentials, but 
rather, after the initial burst from the fast level, molecules in the metastable level exhibit 
an enhanced decay rate, which we have explained to be due to lattice-induced coupling to 
the fast level. This continues for milliseconds and implies that the lattice itself takes ms 
to accommodate to the distortion. This is the slow relaxation, and our ability to use this 
assumption to fit a variety of systems in a large range of temperatures is the evidence for 
this phenomenon. 

Because the distortion of the quasimolecule is asymmetric, a reasonable model is a linear 
diatomic chain in which a particle on one end is subject to a strong push. (This is the Br 
adjacent to the Pb.) Nonlinearity enters because of the large displacements. Moreover, to 
account for the 3-dimensional environment we add a potential (simulating off-chain neigh- 
bors) that tends to return each ion to its normal lattice position. This model displays a 
remarkable property: reverberations in the chain due to the push are confined. Oscillations 
occur only in the 2 to 4 neighbors of the first ion. Except for a small initial pulse traveling 
ballistically, no energy escapes. These assertions are based on numerical integration for 
^ 500 lattice time units. Accurate integration of even the 1-dim classical case for 10 9 units 
(~ms) appears out of reach, so if we claim this confinement underlies slow relaxation, more 
general arguments are needed. 

As we show below, the confinement reflects the creation of a FPU soliton. This structure 
lives almost forever ("almost" refers to Arnold diffusion). We thus go from the problem of fast 
relaxation to that of no relaxation. To explain the ultimate decay of the soliton we note the 
apparent coincidence that the relaxation is on the same scale as the electromagnetic decay 
of the metastable level. This suggests that the electronic degrees of freedom draw energy 
from the soliton. This damping yields patterns of lattice relaxation that follow essentially 
the time dependence used in explaining the experiments. We also found (numerically) that 
even after the quasimolecule decays the soliton remains. 

This explanation of slow relaxation also leads to the question of direct observation. If 
indeed for ~5 ms small bunches of atoms vibrate at a frequency above the acoustical phonon 
branch, what probe could reveal this? 
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Details of the model. Following the UV pulse, and faster than any scale considered here, 
the electronic wave function, if), distorts. The quasimolecule is pushed hard along one axis 
and shrunk along the others. We focus on a ray of atoms along the expansion axis. The first 
atom is subjected to a large force which it transmits to the others. We simplify the influence 
of if) by pretending there is a fictitious zeroth particle displaced by a fixed amount from its 
equilibrium. To allow the use of a 1-dimensional chain we supply each atom on the chain 
with a "neighbor" -force that attracts the atom to its nominal place in the lattice. Note 
that the asymmetric stresses imply lesser deformation for off-axis atoms. The interatomic 
potential is taken to be V(u) = Muq {u 2 + Am 4 ) /2. The Hamiltonian is 

H = £ {H + §| + V(Q n - ftO) + V(Q n - q n ^) + v (V(QI) + V(q 2 n ))} , (1) 

where v is the effective number of neighbors. Q\ is the first atom to the right of the Pb, 
followed by q\, Q2, etc. The Q-particles have mass M, the qs, M/r — but we immediately 
rescale so that M = 1 and Uq = 1. For KBr, r « 2. The impact of if) is expressed by setting 
(the non-dynamic) q 7^ 0. 

The system was solved classically by numerical integration. All positions and momenta 
were initially zero, with q providing the driving force. Fig. |l] shows runs with v = and 4 
(A = 1). In the first case there is wave propagation, but with v = 4 the energy is confined. 
A small pulse leaves the system initially, but afterward only the first atoms oscillate, mostly 
just the first two. Fig. ^| shows the positions for v = 4 but A = 0. Evidently there is 
no confinement. These results are insensitive to the exact parameters and persist both for 
longer times and for more particles. They become more dramatic as r moves away from 
unity, consistent with the trend in anomalous decay data 00]. 
[[ Fig. [l] and Fig. |2| belong here. ]] 
We next study the frequency spectrum of the confined vibrations. To connect to solitons 
this will be related to lattice phonons properties. When v = 4, r = 2, iV = 15, the phonon 
spectrum has a gap between the acoustical and optical branches running from 2.45 to 3.46 
(co>o — 1)- In Fig. ^| we plot the intensity of the Fourier transform of qi{t) for a period of 200 
time units for three values of A. The curve peaked below ~ 2.45 is the linear case, A = 0. 
The next two are A = 0.5 and 1. They are dominated by their peaks, which come at 2.86 and 
3.27 respectively — within the gap. (The bump just below 4 is part of the A = 1 spectrum; 
there is a smaller one above 2.86 for the A = 0.5 case.) If q is varied, the peak-locations 
change, consistent with scaling. 
[[Fig. |3| belongs here. ]] 
Besides the extreme spatial localization, another property of our confinement curves — 
perhaps the characteristic soliton feature in the FPU study — is the failure of the energy to 
disperse among the normal modes. In Fig. [| we show the energy dispersal for various times. 
Unlike FPU it is never concentrated in one mode, because our initial conditions already 
have dispersion. Nevertheless, it is clear that while energy can shift from mode to mode (cf. 
time-25 in the figure), it tends to return (repeatedly). Note that we show only quadratic 
contributions to energy, so the sum displayed need not be conserved. 
[[Fig. |4] belongs here. ]] 
An extensive recent literature on FPU solitons exists |6|-[T0|1, and their presence in alkali 



halide crystals has been suggested f||7j]. Like our confinement mode, these excitations are 
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dominated by a single frequency. In the diatomic chain is studied, and the soliton 
frequency is found between the acoustical and optical phonon bands. 

The ability of our systems to support solitons can be attributed to several features. 
1.) The distortion (a la Jahn- Teller) makes the system closer to 1-d than 3-d, bearing in 
mind J7| the lesser tendency to form solitons in 3-d. 2.) Considerable enhancement of all 
soliton effects develops as |r — 1| f. 3.) From numerical work we find that the boundary 
condition at the impurity site, based on our physical picture, enhances localization. 

Decay of the soliton. Given the confinement -soliton connection, why does the crystal 
relax at all? It is a common misapprehension that the FPU soliton lives forever; it does 
not. In m it is shown that FPU dynamics is equivalent to Henon-Heiles chaos. For our 
parameters, however, as for FPU, the excitation is extremely long-lived, disappearing only 
because of Arnold diffusion. However, even if for KBr:Pb 2+ this time scale would be ms, its 
sensitivity to parameter changes imply that there is no reason to expect the same for other 
doped alkali halides with ms anomalous decay. One should also realize that the soliton 
is robust. Defects may shift its properties but not the fact of its existence. Moreover, 
anomalous decay has been seen with various levels of crystal defects. For these reasons we 
turn to other degrees of freedom, and in view of the apparent coincidence of time scales, it 
is the electronic coupling that we argue provides decay. 

The decay of the metastable A lu level (ip u ) to the A l9 ground state (ip g ) is forbidden by 
electromagnetic selection rules. The way the decay does occur is by coupling %p u to vibronic 
modes, which in turn couple to ip g . We invoke almost the same matrix elements, but in 
a different order: from the solitons to ip u to other lattice vibrations. In one case one has 
(^j|iy|vib.)(vib.|if|^ 9 ), in the other (vib.\H\ip u ) (ip u \H\vib.) . Density of states factors may 
differ. Nevertheless, it is plausible that the rates are comparable, and this is our assumption. 

To incorporate this in our classical model we view the coupling as a damping: energy in 
the solitons is converted to undifferentiated vibrational energy via electronic coupling. This 
affects primarily the Br in the quasimolecule, so that classically we included — 7Q1 (with the 
physical 7 -1 ~ ms). In Fig. |5] we show a smoothing (average) of the oscillations of Q\ over 
an extended time period with damping. This is also superimposed on a curve 1 — exp(— Ft), 
with 7 « T. The reason for this comparison is that 1 — exp(— Ft) is one of the relaxation 
forms adopted in || in the slow-relaxation explanation of anomalous decay. We thus recover 
the desired form for relaxation [|11|] . 

[[Fig. |H belongs here. ]] 
Remarks. 1.) We expect the quantum version our theory to give similar results. In Q 
is a comment on quantum treatments, and we remark that the quantization of the soliton 
can be accomplished by path integral WKB methods \T2\. This is because the phase space 
trajectory of the soliton is essentially a torus. Weak phonon coupling provides damping, as 
in and a lifetime. 2.) For nonzero temperature other mechanisms of soliton decay 
enter and indeed our data fits H indicate increased 'T" with temperature. Mechanisms 
such as those studied in |14| should play a role, although as observed there, the effect may 
be sensitive to details of the system. 

Detection of solitons. Our explanation of slow relaxation implies that for periods ~5 ms 
subsequent to a UV pulse, small bunches of atoms (perhaps 5 on a side) oscillate at fre- 
quencies above the acoustical phonon band. For KBr this is at roughly v a = 3 x 10 12 s~ 1 
|I5fl. A direct way to see these would be by off- resonant Raman scattering, to avoid the 
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difficulties of strong luminescence. The crystal would simultaneously be illuminated with a 
Raman laser and, for soliton creation, UV light (e.g., a hydrogen lamp). One would then 
look for new peaks, in addition to the characteristic peaks of (PbBr) 6 and the second order 
Raman spectrum of KBr. At liquid helium temperatures we expect Raman lines above the 
acoustical branch and therefore any new peak is expected to be near the laser excitation line. 
When the UV light is switched off (terminating A-band excitation of Pb 2+ ) these Raman 
lines should disappear. We also anticipate requiring a higher Pb concentration than is now 
used in decay experiments. 

Finally we observe that the soliton frequency corresponds to an energy, Tiuj ~ lOmeV. 
The energy dumped into vibrational modes is much larger (perhaps as much as 650 meV), 
so that the semiclassical approach to the study of the soliton is reasonable. 
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FIGURES 
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FIG. 1. Position vs. time. The unconfined history has v = 0, no force from off-axis atoms. 
There is wave propagation. The reflection is an artifact of the boundary. The case v = 4 gives 
confinement (for which the boundary is irrelevant). The lowest curve is Qi(t); above it is qi(t), 
etc. Common parameters are A = 1, qo = 1, N = 15 (30 atoms), r = 2. 
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FIG. 2. Position vs. time for a system that with nonlinearity would show confinement. Param- 
eters as in Fig. [j] {y = 4), but A = 0. 
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FIG. 3. Fourier transform (intensity) of Qi{t) for various levels of nonlinearity. The dashed 
curve is A = 0.5. The two solid-line curves do not overlap, with amplitude to the left of 2.45 
corresponding to A = and to the right of 3, to A = 1. For A = 0, the peak is just within the 
acoustical band. With A = 0.5 (dashed line) almost all energy is in the peak, which now lies well 
within the gap. Finally, with A = 1 (solid line, again) the peak has almost reached the bottom of 
the optical branch and a second small peak beyond the optical spectrum is more pronounced. 
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FIG. 4. Energy per normal mode at various times. Parameters as in Fig. [p, v = 4. 



8 




time (delayed) -> 400 500 



FIG. 5. The function Qi(t), partially smoothed, calculated with damping (7). Also shown is 
the curve 1 — exp(— Ft) for T ~ 7 ~ 0.01. Qi(t) is from the confinement run of Fig. [I]. 
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